Plane fronted parallel waves in a warm two - component plasma 
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A general system of equations is derived, using the 1+3 orthonormal tetrad formalism, describing 
the influence of a plane - fronted - parallel gravitational wave on a warm relativistic two - component 
plasma. We focus our attention on phenomena that are induced by terms that are higher order in the 
gravitational wave amplitude. In particular, it is shown that parametric excitations of ion - acoustic 
waves takes place, due to these higher order gravitational non-linearities. The implications of the 
results are discussed. 
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^P' I. INTRODUCTION 

. There have been many theoretical investigations on the scattering of electromagnetic waves off gravitational fields, 
using linearized gravitational wave theory (see e.g. (]]), which is the relevant regime for gravitational wave detectors, 
or, in general, for distances far away from the gravitational wave source. Much of previous research has been directed 
at studying effects on vacuum electromagnetic fields, but there has also been some work where the effects of plasmas 
have been taken into account (see e.g. [§P||] and references therein ). 

In a recent paper ||] we have taken another approach, starting with an exact gravitational wave solution, but 
focusing on a weak amplitude (but still non- linear) approximation, and studying the effects induced by gravitational 
waves in a plasma. As was shown in such interactions give rise to qualitatively new phenomena that are absent 
in linearized theory. However, the previous work was limited to cold plasmas, and thus did not allow interaction 
with low - frequency longitudinal waves (although for example Ref. [||] treats the generation of low - frequency plasma 
modes by linear gravitational waves). As will be shown, such waves are more likely to be excited by gravitational 
waves,[] due to the frequency matching condition. 

We start by generalizing previous equations for a cold two - component plasma to include thermal effects, using 
the 1+3 orthonormal frame formalism. Within this mathematical framework, the governing plasma equations can 
be written in a simple 3 - dimensional form. They consist of Maxwell's equations with additional charge and current 
^P" densities characterizing the gravitational effects and a set of fluid equations for a warm plasma. To facilitate the 
J> , analysis of the non -linear interaction between a plasma and a gravitational wave, we make use of the plane fronted 
parallel (pp) wave solution of Einstein's field equations. We show that for parallel propagation, excitation of ion- 
acoustic waves can only occur if effects second order in the gravitational wave amplitude are included. The growth rate 
of the second order instability is determined, and the threshold value for excitation is estimated. Some applications 
to astrophysics and cosmology are discussed and our results are summarized in the last section of the paper. 
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II. PRELIMINARIES 



The energy - momentum tensor for each particle species is assumed to take the form of a perfect fluid, T ab = 
(p+p)V a V b +pg ab , where [i is the energy density and p the pressure of each fluid. Splitting the energy - momentum 
tensor in time- and spacelike parts, using an orthonormal frame {e a ,a = 0, ...,4}, the particle and momentum 
conservation equations take the form |§] 
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e ( 7 n) + V-( 7 m>) = -yn (T a 0a + T a QO v a + T a 0a v p ) , (la) 
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generalizing the results presented in || to the case of a warm plasma. Here r a b c are the Ricci rotation coefficients with 
respect to the ONF {e a }. We have introduced the notation E = (E a ) = (E 1 ^ 2 ^ 3 ) etc., V = (e 1 ,e 2 , e 3 ). The fluid 
is assumed to have a four - velocity (7, jv), relative to the orthonormal frame {e a , a = 0, ...,4}, where 7 = (1 — v 2 ) -1 / 2 . 
Furthermore, p = qjn is the charge density. 

We follow the approach presented in Refs. 0,0] for handling gravitational effects in Maxwell's equations. With this, 
these equations was presented in their generality in J|,|]]. In order to address the issue of how strong gravitational 
radiation may be involved in generation of EM waves, we look at the pp- waves (for a discussion of this solution, 
see 0), in the special case of a linearly polarized plane wave the spacetime metric takes the following form 

ds 2 = -di 2 + a{uf dx 2 + b(u) 2 dy 2 + dz 2 , (2) 

where u = z — t, and a and b satisfy ab uu + a uu b = 0, and the subscript u denotes a derivative with respect to retarded 
time. Note that we have chosen a vacuum geometry, i.e. we have omitted the influence of the plasma on the metric. 
In order to make physical interpretations simple, we introduce the canonical Lorentz frame 

e = d t , ei = a~ 1 d x , e 2 = b^dy , e 3 = d z . (3) 

With this choice, the gravitational effects in Maxwell's equations are 

p E = -(In ab) u E 3 , (4a) 

p B = -{\nab) u B 3 , (4b) 

j E = -(InbUE 1 - B 2 ) ex - [\na) u {E 2 + B 1 ) e 2 - (lnab) u E 3 e 3 , (4c) 

j B = -{\nb) u {E 2 + B 1 ) ei + (Ina^E 1 - B 2 ) e 2 - (lnab) u B 3 e 3 , (4d) 

which enters Maxwell's equations as 

V • E = p E + p ch , (5a) 

V-B = p B , (5b) 

E-V X B = -j B -j , (5c) 

B + V X E = -j B , (5d) 

where j = J2 P s <ll nv ls the current density, and the overdot stands for time derivative. 

With respect to the frame choice (||), the conservation equations (Q) lead to the following fluid equations: 

d 

— ( 7 n) + V-(jnv) = 7 n(lnafe) u (l - v\\) , (6a) 
d . „\ , _!_ . fd 



ip+P) ( -Q t +v • VJ 7 v+7 Vp + 7« I— +v • Vjp 

= p ch (E + v X B) + 7(/i + p) [Qna) u vxex + (ln6)„w 2 e 2 ] (1 — Uy) + -f(p + p) [{lna) u v 2 + (lnb) u v$\ e 3 , (6b) 

where v\\ = V3 is the velocity parallel to the gravitational wave propagation direction. These equations should be 
satisfied for each particle species. 



III. LONGITUDINAL WAVE EXCITATIONS 

The terms in equations (§)-(||) which generate effects in the direction of propagation of the gravitational wave, i.e. 
the longitudinal terms, are second order in the gravitational wave amplitude. As was previously shown in j|, these 
second order terms can give rise to qualitatively new phenomena compared to the linear regime. However, for the 
process described in the gravitational wave frequency must be equal to the local plasma frequency, which puts 
rather severe constraints on the possible sources for the gravitational waves. Here we will focus on the excitation 
of longitudinal waves with much lower frequencies, specifically ion -acoustic waves. For simplicity, the equation of 



2 



state is assumed to be an non - relativistic isothermal^ i.e. p — ksTn, ks being the Boltzmann constant and T the 
temperature. 

As a consequence of the low temperature, the inertial mass term fi +p « mn and it is this approximation that will 
be used from now on. Furthermore, we consider the regime ln(a&) <C 1. The background is homogeneous in both 
fluid species, and thus we can take the zeroth order electromagnetic fields to be zero. Neglecting the effect of charge 
separation, we obtain the first order solutions to Eqs. (0) 
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where n is the background number density and v t h = \Jk B T jm the thermal velocity. The above equations hold for 
both particle species so that the density perturbations are in general different for the two fluids J^] 

With the above calculation as a prerequisite, we consider the stability of these solutions, i.e. we make the ansatze 



n(t, z) = n a + n gw (u) + h(t, z) , v(t, z) = v gw (u) + v(t, z) , E(t, z) — E(t, z) 



(8) 



We linearize in n, v and E, and keep terms up to order (lna6) u x(n, v or E). The effect of terms quadratic or higher 
order in (\nab) u is to modify the equilibrium solution (Q) by a factor ~ 1 + (lnet6). Thus we neglect the corresponding 
influence on the wave excitation. 

Using the above perturbation scheme, we can combine equations (|6aJ) and ( |6b| ) to obtain 
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Note that & and S in general are different for different particle species. 
Acting with @i@2 on Eq. (]5a|), we obtain 
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E^Qnab) 



(12) 



where the indices 1 and 2 denotes the two fluid species, and w p i,2 = [<?i 2 n o/ m i,2] 1 ^ 2 is the plasma frequency of each 
particle species. Equation (12) holds for any two -fluid plasma with low temperature in both species, and can describe 
excitation of both high-frequency and low - frequency waves. A process involving high-frequency waves in a cold, 
one - component plasma, was investigated in j|. However, in most cases the gravitational wave frequency is lower than 
the electron plasma frequency. Thus, from now on we will focus on three -wave excitation of low - frequency modes 
(i.e. all frequencies ^ w P ,c) in an electron -ion plasma. Note that the gravitational wave acts as our pump wave. We 
are interested in the case of a periodic source of our pump wave, and following we obtain (to second order in an 
expansion in the amplitude h) 



(Inab) = —h 2 cxp[2iuj gw (z — t)] + c.c. , 



(13) 



where c.c. denotes the complex conjugate. 



2 A more general temperature - to - density profile can easily be incorporated into the calculations. 

3 In principle the electric field due to the possible charge separation should be included in Eq. (f?|). However, since this effect 
is proportional to both (hi ab) u and «th, it will not influence our main result (|17|). 
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The electric field perturbation is assumed to have the form 

E = E + (t) cxp[i(k + z - u + t)] + E-(t) exp[i(fc_z - w_t)] + c.c. , (14) 

and the analogous expressions is assumend to hold for n and v. The time variations of the amplitudes are slow, as 
compared to uo± . The wave numbers and frequencies satisfy the matching conditions 

fc+ + k- = 2uj gw , uj + + oj_ = 2w gw . (15) 

However, for low - frequency waves, the phase velocities are generally much smaller than unity, which means that 
equation (nTHl) can be approximated by 



k + = —k- — k , oj + = oj- = uj gw . (16) 

To be consistent, we should simplify the right hand side of Eq. ( |l^ ) as far as possible, using Vth -C 1, which, 
for example, implies S = —nd(lnab)/dt. Inserting the ansatze for h, v and E in Eq. (|l2|), eliminating the density 
variations using Eq. , and applying the approximate matching conditions ( [l6| ) , results in a growth of the amplitudes 
h, v,E oc exp(rt), with the growth rate 



r = 2i dD ^ k) Y\' 

duj sw 
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^gw 2 \ " 22 



-1 , ,2 



where the sum is over particle species and 

£>(u,,fc) = l-E ,„ 2 _i 2 „2 , (18) 



p.s. 



and the index i refers to the ionsQ The scale length 1/fc of the excited modes, obtained from D(iv gw , k) = 0, satisfies 

^ < \ < , (19) 

where c s = ; + (m c /mi)v^ h g] 1 / 2 is the ion - acoustic velocity. 

By including some mechanism of wave damping, a threshold value h t hr for excitation may be determined. The 
most suitable regime for excitation is T c > 7] in which case ion Landau damping is not effective. Then electron -ion 
collisions is the main disspative mechanism, and 



Hhr 



\J ^c-i/wgw , (20) 



which is similar to the results obtained in [Q. Here v c -i is the electron -ion collision frequency. 

Without working out a detailed theory for the saturation mechanism, we note that the wave growth will stop at 
a level v < v t h,i- This is because the ion -acoustic waves become strongly nonlinear at this level, implying efficent 
energy transport away from the originally excited modes. 

IV. SUMMARY AND DISCUSSION 

In the present paper we have generalized previous equations for a cold plasma in the presence of a strong gravitational 
wave by including thermal effects. As is well known, thermal effects are important for low frequency plasma 
phenomena. Since, typically, the time - scales for gravitational waves are long compared to the time - scales of a plasma, 



4 The dispersion function D(uj,k) can be further reduced using lo -C <^ p ,c, e denoting the electrons. Applied to the expression 
for T, we obtain the last approximate equality in Eq. (|17|). 
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this generalization is important from an applicational point of view. The derived equations provide a framework for 
investigating strong gravitational pulse effects in warm multi - component plasmas. It was shown that the equations 
indeed admit generation of ion -acoustic modes, which are not present in the linearized theory (cf. Ref. [||). 

Since our effect is of second order in the gravitational wave amplitude, possible astrophysical applications are most 
likely to be found close to extreme events, such as binary mergers. The expression for the growth rate is almost 
identical to that of Ref. [Q], where two plasmon decay of a pp-wave was considered. However, a fact that make 
the present process more relevant for astrophysical applications is that it in principle can occur for arbitrarily low 
gravitational wave frequency w gw , be. the plasma frequency is allowed to be much larger than u> gw . 

Note that Eq. (f?| implies that the gravitational wave directly induces charge separation, provided that the thermal 
velocities of the particle species are different (which is in general true for ion - electron plasma) . This charge separation 
could be important in the early universe, since it could lead to a weak statistically homogeneous and isotropic electric 
field. Provided the conductivity of the cosmological plasma remains sufficently low, conditions which exist during 
inflation and the subsequent period of reheating, this primordial electric field could survive immediate dissipation and 
could trigger a period of cosmological magnetogenesis through its interaction with linear gravitational waves J|,[|. 
In this way large scale cosmological magnetic field could be generated via physicical processes inherent to plasmas 
and the geometrical nature of spacetime, rather than invoking field theoretical arguments such as the breaking of 
conformal invariance of electromagnetism. 
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